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Landau levels and states of electrons in a magnetic field are fundamental quantum entities 
underlying the quantum Hall and related effects in condensed matter physics. However, the 
real-space properties and observation of Landau wave functions remain elusive. Here we 
report the real-space observation of Landau states and the internal rotational dynamics of 
free electrons. States with different quantum numbers are produced using nanometre-sized 
electron vortex beams, with a radius chosen to match the waist of the Landau states, in a 
quasi-uniform magnetic field. Scanning the beams along the propagation direction, we 
reconstruct the rotational dynamics of the Landau wave functions with angular frequency 
~100GHz. We observe that Landau modes with different azimuthal quantum numbers 
belong to three classes, which are characterized by rotations with zero, Larmor and cyclotron 
frequencies, respectively. This is in sharp contrast to the uniform cyclotron rotation of 
classical electrons, and in perfect agreement with recent theoretical predictions. 
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Classical electrons in a uniform magnetic field propagate 
freely along the field and form confined circular orbits in 
the plane perpendicular to the field. The angular velocity 
of such orbiting is constant and is known as the cyclotron 
frequency. Accordingly, quantum-mechanical eigenstates of a 
scalar electron in a uniform magnetic field are localized in the 
transverse plane and are characterized by two quantum numbers. 
Excluding the longitudinal motion of the electron (as, for 
example, in two-dimensional (2D) condensed- matter systems), 
this leads to the quantization of the energy levels, which are 
degenerate and are characterized by a single quantum number. 
Quantum electron states and their corresponding energy levels in 
a magnetic field were described by Fock , Landau^ and Darwin^ 
in the early days of quantum theory, and are commonly referred 
to as Landau states and Landau levels. 

Landau eigenstates play a key role in various solid-state 
phenomena, such as the diamagnetism of metals, as well as 
quantum Hall, Shubnikov-De Haas and De Haas-van Alphen 
effects^"^. Landau energy levels reveal themselves in quantum- 
Hall conductance plateaus^; they are measured spectroscopically^; 
and recently they attracted enormous attention in relation 
to graphene systems^"^^. However, Landau levels are highly 
degenerate and do not provide information about the actual state 
and spatial distribution of the electron. Moreover, the drift of the 
states in an external random potential blurs the picture in 
condensed-matter systems and it is impossible to observe the 
fast rotational dynamics of electrons in such systems. Although 
considerable progress was achieved recently in Fourier analysis of 
Landau modes (based on a single radial quantum number), 
their real-space properties remain elusive. Thus, the observation 
of spatially resolved Landau eigenstates and their internal 
dynamics remains a challenging problem. 

Landau states can appear not only in condensed-matter 
systems but also for free electrons in a uniform magnetic field. 
Recently, we argued that, allowing free propagation along the 
magnetic field, the Landau states represent non- diffracting 



versions of the so-called electron vortex beams . 
Furthermore, both the radial and azimuthal (vortex) quantum 
numbers of Landau modes crucially determine their properties 
and evolution. Electron vortex beams were predicted^^ and 
recently generated using transmission electron microscopes 
(TEM)i^-^4, and they promise applications in various areas of 
both fundamental and applied physics^"*"^"* (for a review, see 
ref. 35). Here we report the real- space observation of individual 
Landau eigenstates, which are formed by free-electron vortex 
beams in a uniform magnetic field inside a TEM. We measure the 
fast rotational dynamics of electrons within different states and 
reveal their unusual non -classical behaviour. Instead of cyclotron 
orbiting, we observe that Landau electrons rotate with three 
different angular velocities, determined by the vortex quantum 
number. 

Results 

Rotational dynamics of electrons in quantum Landau states. 

The Landau states of an electron in a z- directed homogeneous 
magnetic field B>0 and uniform gauge with azimuthal vector 
potential = Br/ 2 are described by cylindrical vortex wave 
functions 

Here (r, (p^ z) are cylindrical coordinates, m = 0, ±1, ± 2, ... and 
n = 0, 1, 2, ... are th e azimut hal and radial quantum numbers, 
respectively, wb — \/2h /\e\B is the magnetic length parameter, 
are the generalized Laguerre polynomials and kz is the wave 
number of the free longitudinal electron motion. The transverse 
energy of the electron is quantized according to Landau 
levelsi'^'16 

£^ = :5Q(2iV+l), 

where iV = n + (|m| + m)/2 = 0, 1, 2, ... 
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m = 0 
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Figure 1 | Landau states with different azimuthal indices and electron Bohmian trajectories. Grey-scale plots show the transverse probability 
distributions of the states (1) with radial quantunn nunnber n = 0 and azimuthal numbers m = -1, 0, 1 (which correspond to the Landau levels (2) 
with /V = 0, 0, 1). Three-dimensional streamlines of the probability current in equation (3) (that is, Bohmian trajectories) inside these states are shown 
for three different radii r/wB^O.6, 1, 1.4, in the truncated azimuthal range (p\z = oe(0, n), and for the Larmor propagation distance zg(0, v/Q) (v is the 
electron velocity). Different trajectories are marked by different colours for the sake of convenience. The spiralling of the Bohmian trajectories at the 
maximal-intensity radii r^WB^/\m\ (indicated by the orange lines in the bottom planes) corresponds to the m-dependent internal rotational dynamics with 
frequencies <(D> of equation (4). 
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is the principal Landau quantum number and Q = \e\B/2iii^ is the 
Larmor frequency corresponding to the electron charge e— —\e\ 
and mass jHe. Although the classical electron dynamics is 
determined by the cyclotron frequency co^ — lQy it is the 
Larmor frequency that is fundamental in the quantum 
evolution of electrons 

The Landau states (1) resemble z-propagating free-space vortex 
beams^^"^^, which are characterized by the phase factor 
exp(/m(p), circulating m-dependent azimuthal current and 
kinetic orbital angular momentum L^ — tim per electron. 
However, the kinetic orbital angular momentum of the Landau 
states differs significantly from that of free-space beams and is 
determined by the principal quantum number: L^ — 
E_L/Q = h{2N+\)>0 (ref 16). It is this positive angular 
momentum and the corresponding negative magnetic moment 
Mz— {e/2fi^)Lz<0 that are responsible for the diamagnetism of 
electrons predicted by Landau^'^. The nontrivial angular momentum 
of the Landau states appears because the gauge-invariant probability 
current is modified by the presence of the vector potential A 
(ref 37): j = /i^ ^ Im(i/^* Vi/^) - eA . According to this, 
the current in the Landau states (1) becomes 



r \ w 



(3) 



where (p and z are the unit vectors of the corresponding coordinate 
axes. Here the m-dependent azimuthal term originates from the 
free-space vortex current, whereas the second azimuthal term 
describes the contribution from the vector potential A^p — Br/2. 

Equation (3) describes the spiralling of the Landau electron 
about the magnetic-field direction in the sense of Bohmian 



trajectories, that is, streamlines of the probability current^^"'*^ The 
expectation value of the electron's angular velocity, «(r) = v^(r)/r 
(where v = j/|i/^p is the local Bohmian velocity), can be obtained 
from equations (1) and (3): 



fO for m<0 
Q for m = 0 , 
2Q for m>0 



(4) 



where we used <r ^) = /|iAm,n| ^ ^dr j rdr — 

w^^/\m\y for m/0. Equation (4) represents a highly surprising 
result. It shows that the rotation of electrons in a magnetic field in 
the quantum picture is drastically different from the uniform 
classical orbiting. Instead of rotation with a single cyclotron 
frequency cl>c = 2Q, states with negative, zero and positive 
azimuthal indices rotate with zero, Larmor and cyclotron rates, 
respectively. This result is independent of the radial index n and is 
also valid for any superposition of modes with different n and the 
same m as long as <r"^> = w^^/|m|. The probability- density 
distributions and internal Bohmian trajectories in the Landau 
states (1) with n = 0 and m= — 1, 0, 1 (which correspond to the 
degenerate Landau levels (2) with N= 0, 0, 1) are shown in Fig. 1. 
One can see that in the m = 0 mode the trajectory rotation is 
uniform and corresponds to the Larmor frequency Q. At the 
same time, the rotations inside the m/0 modes depend on the 
radius r and coincide with the averaged values < > (equation (4)) 
at the maximal-intensity radii r = \/|m| . As was recently 
demonstrated for photons'*^''*^ Bohmian trajectories can be 
measured experimentally using statistical averaging over many 
identical single-particle events without interparticle interactions. 




Figure 2 | Schematics of the experiment and the beam parameters, (a) A holographic fork mask generates a row of vortex beams with different 
azimuthal indices m = ..., -5, -3, -1, 0, 1, 3, 5, ... (refs 19,20). These beams are focused by a magnetic lens and are studied in the region of maximal 
quasi-uniform magnetic field. The focal plane is shifted few Rayleigh ranges below the observation plane z = 0 to reduce the Gouy-phase rotation^^'^^. 
A knife-edge stop is placed at Zk<0, where it blocks half of each of the beams. Varying the position z^ of the knife edge, we observe spatial rotational 
dynamics of the cut beams propagating to the observation plane (see Fig. 1). (b) Intensity distribution in the \m\ =1 beams. The radius w of the 
focused beams varies slowly with z. In the highlighted range zg( — 80, —30) |im, the beam radius approaches the magnetic radius, w(z)^Wb, and the 
beams acquire the Landau-state properties (see Figs 3 and 4). 



NATURE COMMUNICATIONS | 5:4586 | DOI: 10.1038/ncomms5586 | www.nature.com/naturecommunications 

© 2014 Macmillan Publishers Limited. All rights reserved. 



3 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms5586 



m = -5 



m = -^ m = 0 m=1 



m = 3 



171 = 5 



-30 |im 



Zk=-50 iim 



Zk=-80 iim 















h (fs) 

-400 -300 -200 -100 0 




-60 -40 
Zk(|im) 

Figure 3 | Experimental images of the cut vortex beams and their m-dependent rotations with the propagation distance, (a) Experimental images of the 
row of the cut vortex beams with different m at different positions of the knife edge, (Fig. 2). The scale bar is 50 nm. The opposite inclination of the 
opposite-m states is because of the residual Gouy-phase diffraction effect^^'^^. At the same time, one can see a slow rotation of the m>0 modes with z^, 
while the m<0 states remain motionless. A quantitative analysis of these m-dependent rotations is depicted in b. The azimuthal orientations of the cut 
modes q) (with respect to the extrapolated reference azimuth (/)|^^^q= (/Jq) are plotted versus z^ and the corresponding timescale tk = -^kA (on the top). 
Three lines correspond to the zero, Larmor and cyclotron rotations predicted for the Landau states in equation (4). Error bars include the uncertainty in 
reading, knife-edge roughness and stage positioning. 



The same conditions are realized in electron- optical measurements 
in electron microscopy"*^. 

Experimental measurements. The main goal of our experiment 
is twofold: first, the creation of free- electron Landau states 
(equation (1)) and, second, the observation of their extraordinary 
internal rotational dynamics (equation (4)). To produce the free- 
electron Landau states, we use the fact that they represent non- 
diffracting versions of the free-space Laguerre- Gaussian vortex 
beams *^' . Such electron vortex beams are generated in a TEM 
using a holographic fork mask (a diffraction grating with a 
dislocation) *^'^^. The mask shown in Fig. 2a has a bar/slit ratio of 
1, and it produces beams with different azimuthal indices m = 0, 
±1, ±3, ± 5, ... for different diffraction orders*^. (For other 
bar/ slit ratios, even values of m can also be produced^^.) The 
vortex beams are then focused with a magnetic lens, which has a 
region of a quasi-uniform z-directed strong magnetic field. 
Figure 2 shows a schematic diagram of our experimental set-up 
with the converging vortex beams. We tune the parameters of the 
system, such that in the z-region of a quasi-uniform magnetic 
field the beam radius w{z) comes close to the magnetic radius 
WB\/\ni\ (Fig. 2b). Thus, the vortex beams approximate Landau 
states (1) with n = 0 and different azimuthal indices m. Then, the 
observation of the peculiar rotational dynamics described by 
equation (4) would verify that the beams indeed acquire 



properties of the Landau modes. The more the beam 
radius deviates from WB\/\ni\, the more <r~^) deviates from 
w^^/\m\y changing the electron angular velocity given by 
equation (4). 

To observe the internal rotational dynamics of equation (4) 
and spiralling Bohmian trajectories (Fig. 1) inside the cylindri- 
cally symmetric beams, we borrow a technique successfully 
employed in optics'*^ ""^^ and also recently demonstrated for 
electrons'*^. Namely, we obstruct half of the beam with an opaque 
knife edge stop and trace the spatial rotation of the visible part of 
the beam when the knife edge is moved along the z axis (Fig. 2a). 
Although such truncation of the beam breaks the cylindrical 
symmetry of the initial vortex state, it does not perturb 
significantly the probability currents in the visible part of the 
beam, so that the truncated beam approximately follows the 
internal Bohmian trajectories of the initial cylindrical state. Note 
that in condensed-matter systems the Larmor rotation of 
electrons is fast and cannot be observed as compared with the 
slow motion of the centre of mass in an external potential*^"*^. 
In contrast, for paraxial electrons in a TEM, the 
transverse Larmor dynamics is slow as compared with the 
relativistic longitudinal velocity of electrons: Qw5«v~c. This 
allows mapping of the internal dynamics on the z axis with 
extremely high resolution (in our experiment, the Larmor time 
scale Q ~ * 8 ps corresponded to the propagation distance 
Zl = v/Q^ 1.7 mm (ref. 16)). 
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The experiment was performed according to the above 
approach in a FEI TECNAI F20 TEM at 200 kV acceleration 
voltage (v^OJc). The focusing lens produced the maximal 
longitudinal field B^1.9T, which corresponds to the Larmor 
frequency 120 GHz (using the relativistic mass lie — ypieo) 
and the magnetic radius Wb — 26 nm. Choosing z — Q as the 
observation plane, the region of interest was ze ( — 80, — 30) )im 
(Fig. 2). In this region the magnetic field was uniform up 
to negligible variations ~10~^B in the longitudinal 
component and ~10~^B in the radial component. The knife 
edge was made from a Si crystal. Its position was varied in the 
region of interest (that is, the propagation distance to the 
observation plane was varied by Az^SOjim) to measure the 
rotations of the images of the cut beam (the corresponding 
Larmor angle is S.(p — Az/zl^ 30 mrad, see Fig. 3). Note that the 
focal plane of the beams was set dXz^% jim, that is, few Rayleigh 
ranges below the observation plane, in order to reduce the 
diffractive Gouy-phase rotation of the images"*^'^^ and to 
improve the accuracy of the measurements using sufficiently 
large vortex radii (Fig. 2a). 

Figures 3 and 4 show the results of the experimental 
measurements of the cut vortex modes at different positions 
of the knife edge. In Fig. 3a the images of the modes with 
m = — 5, — 3, — 1, 0, 1, 3, 5 are shown for three values of Z]^. 
Note that the cut edges of the beams with opposite m have 
opposite inclination with respect to the line joining the beams. 
This is the residual Gouy-phase rotation"* '^^ visible at the 
observation plane. At the same time, a slow rotation of the m > 0 
states as a function of Zy, can be detected visually, while the m < 0 
modes do not experience visible rotation. A quantitative analysis 
of the differential rotations of modes with different m is depicted 
in Fig. 3b. One can clearly see three different rates of rotations for 
modes with m<0, m — 0 and m>0, in precise agreement with 
the prediction of equation (4) and in sharp contrast to the 
classical cyclotron orbiting. This confirms that electron vortices 
form Landau states and acquire their peculiar properties in the 
region of interest. We repeated measurements of the mode 
rotations with Zy, with slightly different defocus values and 
holographic masks (including those producing m — ±2 modes). 
For each of such experiments we determined the average value of 
the rotational velocities <«> —v{d(p/dz\^y for different m. The 
results are shown in Fig. 4. One can clearly see that rotational 
frequencies in different measurements fluctuate around the 
theoretical values of equation (4), and their averages over all 
measurements are in very good agreements with the theoretical 
Landau-state behaviour. Note that variations of the rotational 
frequencies in the m/0 modes and robustness of such 
frequencies for the m — 0 mode can be related to the radial 
dependence (independence) of the rotation in the m^O (m = 0) 
Landau states (see Fig. 1). 

Discussion 

To summarize, the extraordinary m- dependent rotational 
dynamics, which is impossible in classical electron propagation, 
reveals the peculiar behaviour of quantum Landau states. 
Although it is commonly believed that the cyclotron rotation of 
electrons underpins Landau states, we have shown that electrons 
can also rotate in the quantum Bohmian picture with either zero 
or Larmor frequencies. (It is worth remarking that the observed 
m-dependent rotation of electrons can be related to the 
Aharonov-Bohm effect^^. Indeed, it is the asymmetry of the 
azimuthal currents in the m>0 and m<0 modes, caused by 
the presence of the vector potential, that is responsible for the 
Aharonov-Bohm phenomenon *^''*^.) We emphasize several 
striking features and the fundamental importance of these 
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Figure 4 | Averaged rotational frequencies for modes with different 
azimuthal indices m. Averaged rotational rates {co} ^v{d(p/dz[^} (such 
as average slopes of the data in Fig. 3b) are shown for different topological 
charges m. Different data points for the same m correspond to different 
series of measurements, and error bars indicate the s.e.m. in each series. 
The solid lines represent frequencies averaged over all measurements, 
while the dashed lines indicate the theoretical values predicted in 
equation (4). The average values and s.e.m. (indicated as shaded bars) 
from all measurements are (co) = (0.09 ± 0.15)Q for m<0, 
<co> =(0.95±0.03)Q for m = 0 and = (2.06 ± 0.14)Q for m>0. 
This verifies the extraordinary rotational dynamics of electrons in Landau 
states, which exhibit zero, Larmor and cyclotron frequencies for the modes 
with m<0, m = 0 and /n>0, respectively. 



results. First, we demonstrated the appearance of quantum 
Landau states within free- electron optics, rather than in a 
condensed-matter system. Second, in contrast to condensed- 
matter experiments and analyses, we separated modes with 
different azimuthal quantum numbers m (some of which 
correspond to the same Landau levels) and showed that this 
index is crucially important for the electron rotation in a 
magnetic field. Third, while the observation of the fast Larmor 
dynamics is currently impossible in condensed-matter systems, 
this rotation becomes measurable for free relativistic electrons. In 
our set-up, this allowed the detection of rotations with frequency 
Q~ 100 GHz, corresponding to the energy difference 
fiQ~100|ieV. Thus, our results provide new insights into the 
fundamental properties of Landau states and pave the way 
towards detailed investigations of their otherwise hidden 
characteristics. 
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